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The present knowledge of the QCD phase diagram based on simulations 
of lattice QCD is summarised. The main questions are whether there is a 
critical point in the QCD phase diagram and whether it is related to a chiral 
phase transition. It is shown that QCD at imaginary chemical potentials 
has a rich phase structure, which can be determined in a controlled way 
without sign problem and which severely constrains the phase structure at 
real chemical potentials. 

PACS numbers: 05.70.Fh,11.15.Ha,12.38Gc 

1. Introduction 

The QCD phase diagram determines the form of nuclear matter under 
different conditions as a function of temperature, T, and chemical potential 
for baryon number, ^b- Based on asymptotic freedom, one expects at least 
three different forms of nuclear matter: hadronic (low hb,T), quark gluon 
plasma (high T) and colour-superconducting (high fiB, low T), as sketched 
in Fig. [1] (left). Whether and where these regions are separated by true 
phase transitions has to be determined by first principle calculations and 
experiments. Since QCD is strongly coupled on scales of nuclear matter, 
a non-perturbative treatment is necessary and Monte Carlo simulations of 
lattice QCD are a natural approach. Unfortunately, the so-called sign prob- 
lem prohibits straightforward simulations at finite baryon density. There 
are several approximate ways to circumvent this problem, all of them valid 
for fi/T < 1 only [U[2] (with quark chemical potential /x = /x#/3): reweight- 
ing, Taylor expansions in [ijT about zero and simulations at imaginary 
chemical potential, where there is no sign problem. The latter can be either 
analytically continued or used as input for the canonical partition function. 
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Fig. 1. Left: Sketch of the QCD crossover line T c (n) vs. the experimental freeze-out 
curve, which has a larger curvature. The remaining transition structure is sugges- 
tive and not yet conclusive in QCD. Right: Comparison of the phase boundary for 
QCD with Nf = 4 staggered quarks on N t = 4 lattices. From [3]. 



As long as ///T< 1, all give quantitatively agreeing results for observables. 
A direct comparison for the calculation of the phase boundary T c {(jl) for 
a theory with four flavours is shown in Fig. [1] (right) [3]. Note that, on 
finite volumes, a transition is always analytic and the phase boundary only 
pseudo-critical. Determining the order of the transition in the thermody- 
namic limit, and hence the existence of a chiral critical point, requires costly 
finite size scaling analyses and is a much harder task. 

The order of the finite temperature phase transition at zero density de- 
pends on the quark masses and is schematically shown in Fig. [2] (left). In 
the limits of zero and infinite quark masses (lower left and upper right cor- 
ners), order parameters corresponding to the breaking of a global symmetry 
can be defined, and for three degenerate quarks one numerically finds first 
order phase transitions at small and large quark masses at some finite tem- 
peratures T c (m). On the other hand, one observes an analytic crossover 
at intermediate quark masses, with second order boundary lines separating 
these regions. Both lines have been shown to belong to the Z(2) universality 
class of the 3d Ising model [H [9] . A convenient observable is the Binder 
cumulant for the order parameter B 4 (X) = {(X - (X)) 4 )/((X - (A)) 2 } 2 , 
with X = ipip or the Polyakov loop. At a second order transition, in the 
thermodynamic limit takes the value 1.604 dictated by the 3d Ising uni- 
versality class. The critical lines bound the quark mass regions featuring a 
chiral or deconfinement phase transition, and are called chiral and deconfine- 
ment critical lines, respectively. The former has been mapped out on Nt = 4 
lattices [10] and puts the physical quark mass configuration in the crossover 
region. The chiral critical line recedes strongly with decreasing lattice spac- 
ing PU 02]: for Nf = 3, on the critical point m^N = 4)/m n (N t = 6) ~ 1.8. 
Thus, in the continuum the physical point is much deeper in the crossover 
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Fig. 2. Left: Schematic phase transition behaviour of Nf =2 + 1 QCD for different 
choices of quark masses at fi — 0. On finer lattices, the chiral critical line moves 
towards smaller quark masses. Middle + Right: Chiral critical surface swept by 
the chiral critical line as /z is turned on. Depending on its curvature, a QCD chiral 
critical point is present or absent. 

region than on coarse lattices, Fig. [2] (left). 

When a chemical potential is switched on, the chiral critical line sweeps 
out a surface, as shown in Fig. [2j According to standard expectations [13], 
for small m u ^, the critical line should continuously shift with fx to larger 
quark masses until it passes through the physical point at /j,e, corresponding 
to the endpoint of the QCD phase diagram. This is depicted in Fig. [2] (mid- 
dle), where the critical point is part of the chiral critical surface. However, it 
is also possible for the chiral critical surface to bend towards smaller quark 
masses, cf. Fig. [2] (right), in which case there would be no chiral critical 
point or phase transition at moderate densities. 

With the currently available methods, there are then two ways of search- 
ing for signals of criticality. One is to fix the quark masses to given, prefer- 
ably physical values and calculate the effect of switching on fi. The other 
is to tune the quark masses to find the base line of the critical surface and 
then follow its behaviour as a function of Li. 



The first step of a calculation consists in determining the phase bound- 
ary. In the last few years such calculations have become realistic for physical 
QCD. Since in this case the zero density transition is a smooth crossover, 
the phase boundary is only pseudo-critical even in the thermodynamic limit, 
and hence observable-dependent. It has a Taylor expansion 



2. The pseudo-critical temperature 




(1) 



The following two calculations with differently improved staggered fermions 
use the chiral condensate as an observable. In H] the curvature was calcu- 
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lated with improved staggered quarks on Nt = 4, 8 lattices with a physical 
strange quark mass. Light quarks were extrapolated to the chiral limit 
assuming 0(4)-scaling, giving k = 0.059(2)(4). In another calculation [5] 
quark masses were fixed to their physical values, so that no extrapolations 
and scaling assumptions had to be made. Simulations on Nt = 6, 8, 10 lat- 
tices were continuum extrapolated to yield k = 0.066(20. The important 
physics observation from these consistent results is that the curvature is only 
about a third of that of the experimentally measured freeze-out curve [6j, 
as shown schematically in Fig. [T] (left). Thus there appears a gap between 
the freeze-out curve and the QCD phase boundary. 

3. Signals for criticality at fixed mass 

Reweighting methods at physical quark masses on Nt = 4 lattices get a 
signal for a critical point at fif, ~ 360 MeV [14] . Quark masses were tuned 
to the mass ratios m^/nip ~ 0.19, m n /rriK ~ 0.27, close to their physical 
values. A Lee- Yang zero analysis was employed, as shown in Fig. [3l For a 
crossover the partition function Z(V, (3, Nt) has zeroes only for unphysical 
complex values of the lattice gauge coupling, whereas for a true phase tran- 
sition the zeroes accumulate and pinch the real axis in the infinite volume 
limit. A caveat of this calculation is the observation that the critical point 
is found in the immediate neighbourhood of the onset of pion condensation 
in the phase quenched theory, which is where the sign problem becomes 
maximal [15J. Therefore, it would be good to have a confirmation with an 
independent method. 

One may also search for a critical point using the Taylor expansion. In 
this case a true phase transitions will be signalled by an emerging non- 
analyticity, or a finite radius of convergence for the pressure series about 
/i = 0, as the volume is increased, to be identified with the critical point, 

P °° ( A 4 A 2n f^E 

F = £c 2n (T)(-J , ^=Jimp n , Pn = 

n=0 

Theorems ensure that if the limit exists and asymptotically all coefficients 
of the series are positive, then there is a singularity on the real axis, which 
would represent a critical point in the (//, T)-plane. The current best at- 
tempt is based on four consecutive coefficients, i.e. knowledge of the pressure 
to eighth order, and a critical endpoint for the Nf = 2 theory was reported 
in [16]. There are also difficulties in this approach. Firstly, there are dif- 
ferent definitions for the radius of convergence, which are only unique in 
the asymptotic limit, but differ quantitatively at finite order. Furthermore, 
estimates at a given order are neither upper nor lower bounds on an actual 
radius of convergence. Finally, finite estimates obtained at finite orders are 
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Fig. 3. Left: Imaginary part of the Lee- Yang zero closest to the real axis |14) . 
Middle: Estimates for the radius of convergence from different observables and at 
different order for Nf = 2 + l,Nt = 4. Dashed lines are hadron resonance gas 
values [17] . Right: First order coexistence region in baryon number density from 
the canonical ensemble on small 6 3 x 4 lattices, ~ 700 MeV [IS] . 



a necessary, but not a sufficient condition for the existence of a critical point. 
For example, one also obtains finite estimates from the Taylor coefficients 
of the hadron resonance gas model, which does not feature a non-analytic 
phase transition. This is illustrated in Fig. [3] (middle) [17j . 

Another signal for a critical point [18j is based on simulations using the 
canonical ensemble of quark number n, which is the Fourier transform of 
the grand canonical partition function evaluated at imaginary ^u, 

Z c (V,T,n) = e-^Ziy^in)^. (3) 

In this case the sign problem is deferred to the Fourier transform and re- 
stricts the total quark number and hence the accessible volumes severely. 
The simulations in [18] are with Wilson fermions run on 6 3 x 4 lattices with 
rather heavy pions, ~ 700 — 800 MeV. A first order coexistence region 
is detected which merges at a critical point, whose baryon density can be 
converted to [J-b/T ~ 2.6, Fig. [3] (right). In this case there is no reweighting 
and no Taylor expansion involved. However, the lattice is coarse, its volume 
very small and the pion mass very far from physical, so again we cannot yet 
conclude for physical QCD. Because of the Fourier transform, the compu- 
tational cost for extrapolations to the thermodynamic limit and to physical 
quark masses are growing nearly exponentially. 



4. The chiral critical surface 

Rather than fixing to a definite set of quark masses, we shall now discuss 
the behaviour of the chiral critical surface in Fig. [2j For definiteness, let us 
consider three degenerate quarks, represented by the diagonal in the quark 
mass plane. The critical quark mass corresponding to the boundary point 
of the chiral transition region has an expansion 

TTvT = 1 + 2^ c *> • V 4 
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Fig. 4. Leftt: Finite difference quotient for the curvature of the chiral critical 
surface, Eq. (J3J) , for Nf = 3 on N t — 6. Middle=Right: The width of the crossover 
region should shrink to zero when approaching a critical point. Instead it is found 
to widen with chemical potential, implying a softening of the crossover [5]. 



Tuning to m c (0), one may evaluate the leading coefficients of this expansion. 
In particular, the sign of c\ will tell us which of the scenarios in Fig. [2] is 
realised. The curvature of the critical surface in lattice units is directly 
related to the behaviour of the Binder cumulant via the chain rule, 

dam c dBi f dB^\ 

d(afj>) 2 d(a[i) 2 \dam J 

In order to guard against systematic errors, this derivative has been eval- 
uated in two independent ways. One is to fit the corresponding Taylor 
series of -B4 in powers of fi/T to data generated at imaginary chemical 
potential, the other to compute the derivative directly and without fit- 
ting |10t 119] . Both methods of calculation give mutually compatible re- 
sults. After continuum conversion one finds for Nf = 3 on Nt = 4 that 
ci = — 3.3(3), C2 = —47(20) [19J. The same behaviour is found for non- 
degenerate quark masses. Tuning the strange quark mass to its physical 
value, we calculated mf d (fi) with c% = -39(8) and c 2 < {20]. Hence, 
on coarse Nt = 4 lattices, the region of chiral phase transitions shrinks, 
i.e. the phase transition weakens at least initially with a real chemical po- 
tential, and there is no chiral critical point for {is ^500 MeV, as in Fig. [2] 
(right). This statement appears robust when a finer lattice is considered. 
As discussed in the zero density section, on Nt = 6 the baseline of the chiral 
critical surface moves closer to the zero mass origin, whereas its curvature 
remains negative, Fig. [J] (left). 

Indeed, the same observation can also be made at fixed masses. At 
zero density, the QCD transition is just an analytic crossover, for which T c 
depends on the observable. In particular, one may use this in order to define 
a width of the crossover region, for details see [5], and study its behaviour 
as a function of chemical potential. If there is a critical point in the phase 
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Fig. 5. Left: phase diagram for imaginary /i. Vertical lines are first order transitions 
between Z(3)-sectors, arrows show the phase of the Polyakov loop. The \i = 
chiral/deconfinement transition continues to imaginary /i, its order depends on Nf 
and the quark masses. Middle: Order of Z{2>) end point for Nf — 2,3 at \x — inT. 
Solid lines are lines of triple points ending in tricritical points, connected by a Z{2) 
critical line. Right: Generalisation of the previous to non-degenerate quark masses. 

diagram, the definition of T c becomes unique and the lines T c (/x) computed 
from different observables must meet at a critical point. In [5] the width 
of the crossover region was evaluated based on the curvature of T c , i.e. its 
leading order Taylor coefficient, and extrapolated to the continuum, with a 
result shown schematically in Fig. 01 The width grows slightly rather than 
getting narrower, i.e. the crossover gets even softer initially as real chemical 
potential is switched on. 

Note that one also observes a weakening of the phase transition with \x 
in the heavy quark case [U] as well as a weakening of the transition with 
isospin chemical potential [21], suggesting that this is a generic feature for 
gauge theories with chemical potential for fermion number. 

5. Critical surfaces at imaginary chemical potential 

Recent Monte Carlo studies at imaginary chemical potential indicate 
that at least for QCD this appears to be the general behaviour, and that the 
sign of the curvatures of the critical surfaces can be understood. Moreover, 
since simulations at imaginary chemical potential have no sign problem, 
they are completely controlled. As we shall see, there are interesting phase 
structures which can then be used to constrain model or effective theory 
descriptions of QCD. 

The QCD partition function is cyclic in imaginary chemical potential, 
Z(/i/T) = Z(fi/T + i2irn/3), due to its Z(3) centre symmetry. This implies 
transitions between adjacent centre sectors, distinguished by the phase of 
the Polyakov loop, at imaginary chemical potentials [l\ = (2n + 1)ttT/3. 
The schematic phase diagram is shown in Fig. Transitions in \i{ between 
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neighbouring sectors are of first order for high T and analytic crossovers 
for low T \2'1\ I23j . as shown in Fig. [5] (left). Correspondingly, for fixed 
fii = fj.%, there are transitions in T between an ordered phase with two-state 
coexistence at high T and a disordered phase at low T. An order parameter 
to distinguish these phases is the shifted phase of the Polyakov loop, L = 
\L\ exp(ic^), 4> = ip — Hi/T. At high temperature it fluctuates about (<fi) = 
±7r/3 on the respective sides of \x\. The thermodynamic limit picks one of 
those states, thus spontaneously breaking the reflection symmetry about [x\. 
At low temperatures <j) fluctuates smoothly between those values, with the 
symmetric ground state {(f)) = 0. Away from fii = fj,f, there is a chiral or 
deconfinement transition line separating high and low temperature regions. 
This line represents the analytic continuation of the chiral or deconfinement 
transition at real [i. Its nature (1st, 2nd order or crossover) depends on 
the number of quark flavours and masses. The junction between the Z{3) 
transition and the chiral/deconfinement transitions at fixed [i = mT/3 was 
studied in [JH [25] for Nf = 2, 3, respectively, with qualitatively similar 
results. For small/large quark masses the first order chiral/deconfinement 
transition connects to the Z(3) transition, and the junction is a triple point. 
For intermediate quark masses they do not connect, and the Z{3) transition 
has a Z{2) end point. This results in the phase diagram for fixed shown 
in Fig. [5J (middle), which is qualitatively the same for Nf = 2, 3. 

The simplest generalisation to non-degenerate quark masses is then 
shown in Fig. [5] (right), where the respective tricritical points for Nf = 2, 3 
are smoothly connected by tricritical lines, both in the heavy and light mass 
regimes. Note that this diagram is the analogue for [i = iirT/3 of the di- 
agram Fig. [2] (left) for fi = 0. The key observation is now that the chiral 
and deconfinement critical surfaces continue to imaginary and terminate 
in the tricritical lines |25j . as shown in Fig. [6] (left). In the case of heavy 
quarks, the critical surface is known over a large range of imaginary and 
real fi within an effective theory [9], the 3d Z(3) Potts model, which is in 
the same universality class as QCD with heavy quarks. For fixed flavour 
content, i.e. a slice through the critical surface, the deconfinement critical 
quark mass follows tricritical scaling [25] ; Fig. (right), 



The shape of the deconfinement critical surface is thus determined by the 
tricritical scaling law, while the sign of its curvature at [i = follows from 
the fact that m tric {ij, = inT/3) < m c (/i = 0). While the chiral critical 
surface is not yet mapped out, one finds in the light quark mass regime 
mtricin = iirT/3) > m c (/i = 0), which thus favours a negative sign for its 
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Fig. 6. Left: The critical surfaces delimiting the chiral and deconfmement transi- 
tions continue to imaginary, where they end in tricritical lines. Right: Critical line 
TO c 2 ) in the 3-state Potts model fitted to Eq. © [3 US- 
curvature. These findings are consistent with a monotonous weakening of 
the chiral and deconfinement transitions as fi 2 gets more positive, which 
we now understand as being induced by the critical structure at imaginary 
chemical potential. 



In summary, signals from reweighting, radius of convergence estimates 
and canonical methods are consistent with a critical point, but their sys- 
tematics does not yet permit definite conclusions for physical QCD in the 
continuum. On the other hand, following the chiral critical surface with 
comparatively controlled systematics tells us that the chiral phase transi- 
tion weakens with moderate ^, thus leading us away from the physical point. 
Possible scenarios are: if we mistrust the systematics of the former meth- 
ods, we would conclude for either no critical point at all, or a critical point 
at larger chemical potential /X£>500 MeV, where current methods break 
down. However, another possibility is that all calculations hold, with a crit- 
ical point at moderate densities that would not belong to the chiral phase 
transition, but to physics unrelated to chiral symmetry breaking. It would 
be interesting to distinguish these possibilities with the help of effective 
theories, which can be gauged against the critical structures at imaginary 
chemical potentials. 



6. Conclusions 
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